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The process of diffusive acceleration of charged particles in shocked plasmas is widely invoked in astrophysics 
to account for the ubiquitous presence of signatures of non-thermal relativistic electrons and ions in the universe. 
A key characteristic of this statistical energization mechanism is the absence of a momentum scale; astrophysical 
systems generally only impose scales at the injection (low energy) and loss (high energy) ends of the particle 
spectrum. The existence of structure in the cosmic ray spectrum (the "knee") at around 3000 TeV has promoted 
contentions that there are at least two origins for cosmic rays, a galactic one supplying those up to the knee, 
and even beyond, and perhaps an extragalactic one that can explain even the ultra-high energy cosmic rays 
(UHECRs) seen at 1-300 EeV. Accounting for the UHECRs with familiar astrophysical sites of acceleration has 
historically proven difficult due to the need to assume high magnetic fields in order to reduce the shortest diffusive 
acceleration timescale, the ion gyroperiod, to meaningful values. Yet active galaxies and gamma-ray bursts remain 
strong and interesting candidate sources for UHECRs, turning the theoretical focus to relativistic shocks. This 
review summarizes properties of diffusive shock acceleration that are salient to the issue of UHECR generation. 
These include spectral indices, acceleration efficencies and timescales, as functions of the shock speed and mean 
field orientation, and also the nature of the field turbulence. The interpretation of these characteristics in the 
context of gamma-ray burst models for the production of UHECRs is also examined. 



1. INTRODUCTION 

The origin and nature of very high energy and 
ultra-high energy cosmic rays (UHECRs) contin- 
ues to pique the interest of the physics and as- 
trophysics communities. This topicahty is driven 
by the lack of resolution of the apparent incom- 
patibility of the data taken by the AGASA and 
Fly's Eye experiments at energies greater than 
around 6 x 10^^ eV. Moreover, both these exper- 
iments and the Yakutsk initiative claim detec- 
tions [37116155] of cosmic ray events above the 
so-called Greisen-Zatsepin-Kuzmin (GZK) cut- 
off, predicted by Greisen 36 and Zatzepin and 
Kuzmin 66 . This observational characteristic 
can place stringent constraints on both the seed 
for production of such energetic particles, and the 
sites in the Universe for their generation. 

The UHECR paradigm is also underpinned by 
the existence of two competing scenarios for their 
creation: (i) the more conventional bottom-up 
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models, where the particles are accelerated from 
much lower energies up to the extreme values in- 
dicated by the air shower array data, and (ii) the 
generally newer and more exotic "top-down" pos- 
sibilities that invoke creation and decay of various 
entities in particle physics that cascade down in 
energy to the GZK domain. Top-down scenarios 
have emerged as competitors to the more tradi- 
tional bottom-up concepts since the latter have 
difficulty achieving energies in the 10^" eV range 
in most known astrophysical objects. 

This review focuses on the properties of dif- 
fusive acceleration at relativistic shocks that are 
relevant to the production of UHECRs in astro- 
physical settings. These characteristics are not 
widely known outside the cosmic ray accelera- 
tion field. The viability of any putative source 
of UHECRs is contingent on its ability to gener- 
ate particles of requisite energies in the available 
time, and with sufficient abundances and appro- 
priate metallicity. Such quantities are sensitive to 
assumed parameters of shocks, such as field obliq- 
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uity, level of turbulence, the nature of scattering, 
and the role of accelerated particles in shaping 
the shock hydrodynamics. Subtle (and even not 
so obscure) implications of acceleration proper- 
ties are often not addressed in bottom-up models 
for UHECRs that are offered for peer scrutiny, 
being left to specialist discussions within the ac- 
celeration community. 

Here, the goal is to highlight some key prop- 
erties of relativistic shock acceleration that are 
pertinent to various UHECR scenarios. To pre- 
pare the way, a few standard results from acceler- 
ation theory at non-relativistic shocks are briefly 
reviewed in Section [21 to facilitate comparison 
with the relativistic domain. This is the most- 
studied aspect of shock acceleration theory, with 
a number of very instructive reviews 23 17 41 50 
in the literature. Moreover, it is the eminently 
testable domain, affording observational diagnos- 
tics via in situ spacecraft measurements 33 10 
of accelerated populations at shocks in the he- 
liosphere, and also numerous particle production 
and radiation models of supernova remnants (see 
[H] for a review), the principal contender for the 
sources of galactic cosmic rays. Section 13 then 
turns to the intricacies of relativistic shocks, fo- 
cusing on several characteristics that distinguish 
them from their non-relativistic counterparts. A 
discussion of the impact of these properties on the 
gamma-ray burst (GRB) paradigm for UHECR 
production is offered in Section 0] 

2. NON-RELATIVISTIC SHOCK AC- 
CELERATION THEORY 

The theory of diffusive acceleration of particles 
at non-relativistic shocks (those with upstream 
fluid speeds ui ^ c in the shock rest frame) 
has been thoroughly investigated by both ana- 
lytic and simulation techniques. The process, also 
know as Fermi acceleration, proceeds by particles 
being transported back and forth across a shock 
between colliding fluid flows. The kinematics of 
momentum deflections then yields a net energy 
gain, and the process can continue for a large 
number of shock crossings to generate high en- 
ergy cosmic rays before they are lost due to con- 
vection, upstream escape, or even cooling. Much 



is understood about this process, though the key 
outstanding issue is how electrons and ions are 
heated in the shock layer, and subsequently in- 
jected into the diffusive acceleration process. It 
is instructive to review three basic characteristics 
of Ml ^ c shocks to set the scene for the rela- 
tivistic considerations of Section |31 

2.1. Canonical Index 

Non-relativistic shocks generate particles with 
a power-law distribution in momentum. This is a 
consequence of high energy particles (those with 
speeds v ^ ui) attaining isotropy in all perti- 
nent reference frames, so that the so-called diffu- 
sion approximation can be applied. At such ener- 
gies, the principal transport equation describing 
the acceleration process, the diffusion-convection 
equation, can be solved analytically for plane 
shocks |19I41| . yielding the well-known result for 
the momentum distribution 



where r = ui/u2 is the shock (velocity) com- 
pression ratio, p is the momentum. Eq. (^) 
is a steady-state, test-particle result. In this 
limit, the spectral index, a , is independent of the 
shock speed, ui , the field obliquity, and any de- 
tails of the scattering process as long as isotropy 
of highly super-thermal particles is maintained. 
The canonical nature of this result is a driving 
force behind invocations of acceleration in astro- 
physics. Note that a high Mach number, non- 
relativistic shock has r 4 and dn/dp oc . 

2.2. The Character of Oblique Shocks 

The result in Eq. ^ exhibits no information 
concerning the normalization of the power-law. 
Such a property is controlled by the rate at which 
particles are injected from thermal energies. This 
injection rate is sensitive to the angle ObhI the 
mean upstream magnetic field makes to the shock 
normal. In the downstream region, particles are 
swept away more efficiently from the shock by the 
convective force of the flow when the obliquity an- 
gle Osni is higher. For non-relativistic shocks of 
high Mach number, when ObhI ^ 30° , the con- 
vection is so effective that injection of thermal 
particles into the acceleration process is entirely 
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suppressed [5]. The critical angle at which accel- 
eration ceases increases to around ObuI ^ 55° 
for hotter, low Mach number shocks. These re- 
sults apply if the diffusion is dominant along the 
magnetic field, i.e. the ratio of spatial 

diffusion coefficients K|| along the field and kj_ 
orthogonal to B is much less than unity. Here 
'^\\,± — '^||,-L''^/3 for scattering mean free paths 
A||^(w), for particle speeds v. A concomitant 
effect of obliquity is the dramatic increase in the 
acceleration rate j39l57j in highly oblique shocks, 
above the ObhI — 0° values [Hi] that can be in- 
ferred from Eq. (O below. This increase can be 
ascribed to the presence of electric fields in the 
shock layer that seed shock-drift acceleration. 

While the rapid reduction of the acceleration 
time in oblique and quasi-perpendicular shocks is 
enticing, the inefficient acceleration is unattrac- 
tive to cosmic ray production models. This trade- 
off is an unavoidable property of oblique shocks 
Restricting to k± <Si k\\ cases is not al- 
ways appropriate. Increasing k± , physically cor- 
responding to stronger field turbulence, tends to 
eliminate laminarity information in the field and 
accordingly renders the shock more like a paral- 
lel, ObiiI ~ 0° one. Consequently, increasing k± 
towards the Bohm limit of k± ~ K|| mutes the in- 
creases in acceleration rates, and returns injection 
efficiencies to levels commensurate with plane- 
parallel shocks Modelers therefore must 
compromise: either opt for quasi-laminar oblique 
shocks that are faster but less efficient acceler- 
ators, or extremely turbulent shocks of arbitrary 
obliquity that generate cosmic rays efficiently but 
at standard rates, i.e. the inverse gyrofrequency 
corresponding to a specific particle energy. 

2.3. Nonlinear Modifications 

The aforementioned features apply to circum- 
stances where the accelerated particles do not 
modify the shock hydrodynamics, i.e. are dynam- 
ically unimportant and act as test particles. Yet, 
when acceleration is extremely efficient, a sizable 
fraction of the total energy budget emerges as 
high energy cosmic rays, an inevitable occurrence 
in an r « 4 shock if the power-law in Eq. Q 
extends to high enough energies. The pressure 
of these particles decelerates the upstream flow, 



which in turn provides feedback to the distribu- 
tion of accelerated ions and electrons, and there- 
fore the fraction of energy going into these par- 
ticles. This nonlinearity has been thoroughly ex- 
plored in the literature (24.23 2 5 31 27 1 5 49l^ 
and is a critical characteristic of efficient shocks. 
The quintessential example is the Earth's bow 
shock immersed in the solar wind; it affords a nice 
data comparison between experiment and accel- 
eration theory 

As a result of the energy conservation that 
regulates the acceleration and the energy appor- 
tionment between thermal ions and high energy 
cosmic rays, there are two distinctive features of 
nonlinear, cosmic ray modified shocks: (i) a dis- 
tribution that deviates from pure power-law na- 
ture, exhibiting a characteristic upward concav- 
ity due to higher energy particles sampling larger 
effective compression ratios, since their diffusive 
mean free paths are longer, and (ii) the thermal 
particles are somewhat cooler ^3] than for test- 
particle shocks, since the subshock is weakened 
and energy removed from the thermal ions and 
electrons. These phenomena can be probed to 
a certain extent by examining isolated, radiat- 
ing systems such as supernova remnants, and at 
present, there are at best modest indications to 
support these theoretical predictions. 

3. DISTINGUISHING PROPERTIES OF 
RELATIVISTIC SHOCKS 

Diffusive acceleration at relativistic shocks is 
far less studied than that for non-relativistic 
flows, yet it is a most applicable process for 
UHECR generation, and may occur in extreme 
objects such as pulsar winds, hot spots in radio 
galaxies, jets in active galactic nuclei and micro- 
quasars, and GRBs. Early work on relativistic 
shocks was mostly analytical in the test-particle 
approximation (e.g., ^5.9 45 38 47 J, although the 
analytical work of |b(JI9j explored nonlinear, cos- 
mic ray modified shocks. Complementary Monte 
Carlo techniques have been employed for rela- 
tivistic shocks by a number of authors, includ- 
ing test-particle analyses by |4(il32j for parallel, 
steady-state shocks, and extensions to include 
oblique magnetic fields by |58I3I13| . 
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A key characteristic that distinguishes rela- 
tivistic shocks from their non-relativistic coun- 
terparts is their inherent anisotropy due to rapid 
convection of particles through and away down- 
stream of the shock, since particle speeds v are 
never much greater than the downstream flow 
speed U2 ~ c/3 . Accordingly, the diffusion ap- 
proximation, the starting point for virtually all 
analytic descriptions of shock acceleration when 
Ml <^ c, cannot be invoked since it requires nearly 
isotropic distribution functions. Hence analytic 
approaches prove more difficult when 71 ^ 1 , 
though advances in special cases such as the limit 
of extremely small angle scattering (pitch angle 
diffusion) are possible |45I44| . Let us explore 
some of the distinctive properties of particle ac- 
celeration at relativistic shocks. 

3.1. Non-Universality of the Spectrum 

The most attractive feature of non-relativistic 
shock acceleration theory is that the distribution 
of accelerated particles is scale-independent, i.e. 
a power-law, as in Eq. with an index a that 
depends only on the velocity compression ratio 
r = ui/u2, i.e. hydrodynamic quantities. This 
elegant result does not carry over to relativistic 
shocks because of their strong plasma anisotropy. 
As a consequence, while power-laws are indeed 
created, the index a becomes a function of the 
flow speed, the field obliquity, and the nature of 
the scattering, all of which intimately control the 
degree of particle anisotropy. 

In the specific case of parallel, ultrarelativis- 
tic shocks, the analytic work of Kirk et al. |H) 
demonstrated that as Fi 00 , the spectral 
index a asymptotically approached a constant, 
a 2.23 , a value realized once Fi ^ 10 . 
This enticing result, which has been confirmed 
by Monte Carlo simulations [131411125] . has been 
referred to sometimes as an indication of the uni- 
versality of the index in relativistic shocks. In this 
subsection, it is illustrated that the asymptotic 
index of 2.23 is indeed not canonical, but rather 
a special case corresponding to compression ratios 
of r = 3 and the particular assumption of small 
scattering (pitch angle diffusion), specifically for 
incremental changes ^scatt in a- particle's momen- 
tum with angle ^scatt ^ 1/ri . 



First, the spectral index of the power-law dis- 
tribution is a declining function of the Lorentz 
factor for a fixed compression ratio, a charac- 
teristic evident in |45l2l44j for the case of pitch 
angle scattering, and a property that extends to 
large angle scattering |32I4| . Faster shocks gen- 
erate flatter distributions if r is held constant, a 
consequence of the increased kinematic energiza- 
tion occurring at relativistic shocks. Note that 
imposing a speciflc equation of state such as the 
Jiittncr-Synge one renders r a function of Fi so 
that this monotonicity property can disappear, 
as evinced in Fig. 2 of Kirk et al. [SI- Table HI 
lists indices a of the dn/dp distribution obtained 
from the Monte Carlo simulation technique of El- 
lison et al. 121]. These results, obtained specifi- 
cally in the limit of pitch angle diffusion, illustrate 
the ffattening as Fi increases. 



Table 1 

Spectral Indices a for Pitch Angle Diffusion at 
Relativistic, Plane-Parallel Shocks. 
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The spectral indices are from the Monte Carlo 
simulation of [22], and the accuracy of their 
determination is typically of the order of 
±0.02. Notes: (a) Here Pi = ui/c and Fi 
are the dimensionless velocity and Lorentz 
factor of the upstream fiow in the shock rest 
frame, respectively. (b) This non-relativistic 
limit approximately reproduces the well-known 
(T = (r -f 2)/(r - 1) resuh in Eq. ((TJ). 

The choice of the canonical compression ratio 
r = 3 is a well-known result for a relativistic, 
purely hydrodynamic shock possessing an ultra- 
relativistic equation of state ^H]- However, one 
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can envisage situations where the magnetic field 
becomes dynamically important. The classic ex- 
ample is the termination shock for the Crab pul- 
sar wind, where Kennel & Coroniti 143' observed 
that strong fields can weaken magnetohydrody- 
namic shocks considerably. In an interesting gen- 
eralization of this, Double et al. [221 recently de- 
termined deviations from r = 3 in ultrarelativis- 
tic shocks, in cases where pressure anisotropy is 
significant, a characteristic that is expected to be 
common in relativistic shocks. Such departures 
can either strengthen or weaken the shock de- 
pending on the nature of the pressure anisotropy, 
which must be a significant function of the shock 
obliquity, i.e., ObiiI • Hence, we anticipate that 
a will be a function ObhI , an issue visited again 
in this Section. 

More novel is the fact that the slope of the non- 
thermal particle distribution depends on the na- 
ture of the scattering, a feature evident in the 
works of Refs. |32I13I4| . The asymptotic, ul- 
trarelativistic index of 2.23 is realized only in 
the mathematical limit of pitch angle diffusion 
(PAD), where the particle momentum is stochas- 
tically deflected on arbitrarily small angular (and 
therefore temporal) scales. In practice, PAD re- 
sults when the scattering angle ^scatt is inferior to 
the Lorentz cone angle l/Fi in the upstream re- 
gion. In such cases, particles diffuse in the region 
upstream of the shock only until their angle to the 
shock normal exceeds around l/Pi . Then they 
are rapidly swept to the downstream side of the 
shock. The energy gain per shock crossing cycle is 
then of the order of a factor of two, simply derived 
from relativistic kinematics |35l4j . The spectrum, 
depicted in Fig.^ is then slightly steeper than the 
result for a strong, non-relativistic shock, due 
to the balance between particle energization and 
loss by convection downstream. 

The results in Fig. ^ are from the Monte Carlo 
simulation of acceleration at relativistic shocks 
developed by Ellison et al. who demon- 

strated that for large angle scattering (LAS, with 
^scatt TT ) the spectrum is highly structured and 
much flatter than E^'^ . Such a case is exhibited 
in the Figure. The structure is kinematic in ori- 
gin, where large angle deflections lead to distribu- 
tion of fractional energy gains between unity and 
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Figure 1. Particle distributions from a paral- 
lel ( ObhI = ) relativistic shock of r = 3 and 
Lorentz factor Fi = 5 , obtained from a Monte 
Carlo simulation |32I4| . Scattering is modeled by 
randomly deflecting particle momenta by an an- 
gle ^scatt within a cone whose axis coincides with 
the momentum prior to scattering. Distributions 
are depicted for three cases, ^scatt < 14° , cor- 
responding to pitch angle diffusion (PAD), large 
angle scattering (LAS: 6'scatt < tt » 1/Fi ), and 
an intermediate case (dotted histogram). 



F^ . Gains like this are kinematically analogous 
to the energization of photons by relativistic elec- 
trons in inverse Compton scattering. Each struc- 
tured bump or spectral segment in Fig. ^ corre- 
sponds to an increment in the number of shock 
crossings, successively from 1 — > 3 — s- 5 — > 7 
etc., as illustrated by Baring T, that eventually 
smooth out to asympotically approach an index 
of a ^ 1.5. Clearly, such highly-structured dis- 
tributions have not been inferred from radiation 
emission in any astrophysical objects. 

An intermediate case is also depicted in Fig.^ 
with ^scatt ~ 4/ri . The spectrum is smooth, like 
the PAD case, but the index is lower than 2.23. 
Astrophysically, there is no reason to exclude such 
cases. Moreover, from the plasma point of view, 
magnetic turbulence could easily be sufficient to 
effect scatterings on the order of these angles, a 
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contention that becomes even more salient for ul- 
trarelativistic shocks with Fi 10 . Note that 
the Fi = 5 results depicted here are entirely rep- 
resentative of the nature of such ultrarelativistic 
cases. Clearly a range of indices can be supported 
when ^scatt is chosen to be of the order of 1/Fi , 
and the scattering corresponds to the transition 
between the PAD and LAS limits. It is antici- 
pated that various astrophysical systems will en- 
compass a range of scattering properties. Accord- 
ingly, the continuous and monotonically decreas- 
ing behavior of a with ^scatt , as indicated in the 
exposition of '29', highlights the non- universality 
of the distribution index in rclativistic shocks. 

Categorizing the scattering as either PAD or 
LAS is a useful division, but is not a complete de- 
scription of diffusive transport in shocks. Another 
characteristic, the diffusion of particles across 
mean field lines, becomes a critical element in the 
discussion of oblique or perpendicular shocks. As 
mentioned above, when the upstream angle 0Bni 
of the field to the shock normal is significant, 
diffusion of particles in the downstream region 
struggle to compete with convective losses, and 
transport back upstream of the shock layer be- 
comes inefficient. In non-relativistic shocks, this 
effect was explored by Baring et al. ,8 , who 
observed that the losses controlled the injection 
efficiency of nonthermal particles, so that when 
©Bnl ^ 30° , thermal particles of speed v ^ ui 
fail to return to the shock after one crossing to the 
downstream side and the Fermi acceleration pro- 
cess is quenched. Accordingly, for ObuI < 30° in 
such regimes, while power-law superthermal dis- 
tributions are exhibited, their normalization is a 
strongly declining function of ©Bni when trans- 
port across field lines is suppressed. 

This phenomenon is manifested in a somewhat 
different manner in rclativistic shocks. When 
Ml '~ c , the V ^ ui criterion for dramatic, if 
not catastrophic, convective losses is satisfied for 
all particle speeds, not just slightly suprathermal 
ones. Hence such losses can be expected to be 
pervasive for all non-thermal energies. Increased 
losses must diminish the nonthermal population, 
and since the loss rate is purely a function of par- 
ticle speed j59l41) . which is effectively pinned at 
V K, c , and U2 , the net effect is to increase the 



spectral index while retaining power-law charac- 
ter. This property is illustrated in Fig. |21 where 
the simulation output was acquired in the absence 
of cross field diffusion (i.e., k± =0). Increasing 
©Bnl results in a rapid rise in <t corresponding 
to a suppression of acceleration. Essentially, for 
©Bnl 25° , acceleration is virtually non-existent 
for Fi/3i ^ 1 . 
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Figure 2. Particle distribution indices a from 
oblique ( ©Bni > ) rclativistic shocks of r = 3 
and different Lorcntz factors Fi , obtained from 
a Monte Carlo simulation |32I4| in the limit of 
pitch angle diffusion. Results are depicted for 
the case of zero diffusive transport perpendicular 
to the mean field, i.e., — Q for the perpen- 
dicular component of the spatial diffusion coef- 
ficient K, . The index is insensitive to ©Bni for 
non-relativistic shocks, but rapidly increases with 
obliquity for rclativistic ones, underlining their 
inherent inefficiency. 



Turbulent plasmas in shock environs generally 
will not admit a = assumption. Strong 
turbulence will drive the system towards the so- 
called B ohm- diffusion limit, where diffusion co- 
efficients are similar parallel and perpendicular 
to the field, i.e. k,±_ ^ K|| , and transport is 
effectively isotropic. Efficient transport across 
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field lines can to a significant extent circumvent 
convective losses, returning the particles to the 
shock from the downstream region, and accord- 
ingly flattening the power-law distribution. It is 
anticipated that transport near the Bohm limit 
would be essential to generate ct 5s 3 , i.e. indices 
meaningful for UHECR acceleration paradigms. 
This is borne out in ,13) and the recent work of 
Ellison & Double E^l, who obtained a « 2.34 
for a Fi = 10 , ObuI = 60° shock in the ex- 
treme case of the Bohm limit. The Monte Carlo 
simulation results of |13l29j also exhibited the ex- 
pected monotonic decrease in cr with the increase 
in . 

3.2. Acceleration Times 

Having explored spectral issues germane to the 
UHECR problem, we now turn our attention to 
the maximum energy issue. This is essentially de- 
termined by the rate at which particles are accel- 
erated, so diffusive acceleration times become the 
focal point. Various authors have researched this 
subject for relativistic shocks 32 12 11 .1.6 52^. In 
particular, [221 found that for large angle scatter- 
ing, the acceleration time for a Fi ^ 5 shock was 
only marginally shorter than that expected from 
classical non-relativistic shock theory. The sim- 
ulations of Bednarz & Ostrowski I 12lll| revealed 
similar modest reductions for both LAS and PAD. 
Observe that in spite of substantial energy gains 
per shock crossing, typically on the order of F^ , 
the particles then spend considerable time diffus- 
ing downstream, a time coupled to their inverse 
gyrofrequency. 

Recently, Baring computed acceleration 
times in the limit of pitch angle diffusion using 
the simulation of _32 , specifically for application 
to jets in blazars. It was found that extrapolation 
of simulations into the relativistic regime revealed 
a hard lower bound on the total acceleration time 
Tacc IS measured in the shock rest frame. The 
time Tacc monotonically decreases (for ultrarela- 
tivistic particles) to this limit as Fi increases to 
infinity, yet proximity is achieved for Fi ^ 10 . 
If Vg represents the energy-dependent gyrofre- 
quency of an ultrarelativistic electron or ion, then 
the velocity dependence of the acceleration times 
in plane-parallel shocks can be approximated (to 



around 1-3% accuracy) by the empirical fit 

(l_ 0.18 1 0.22 ^ 

^""^ I 4 Fi/3i +rf^+ l + Fi/3i 

(2) 

f 

TnR = 7-NR,(/3i = 1) . 

Here Tj^niPi = 1) is the extrapolation of the 
well-known acceleration time formula |34I39) for 
non-relativistic, parallel shocks to fiow speeds c . 
The times are for a velocity compression ratio of 
Ui/u2 = 3 , and the coefficient / describes de- 
tails of the differences in diffusion between the 
upstream and downstream regions, and is of the 
order of unity and independent of Fi . Note that 
when /3i <g; 1 , the familiar non-relativistic re- 
sult emerges: Tacc = Tnh,/ Pi ■ Introducing shock 
obliquity can speed up the acceleration, as in non- 
relativistic cases, but at the price of dramatically 
steepening the distribution. 

The bound arises due to the insensitivity of 
the downstream fiow speed and diffusion in the 
downstream region to the upstream Fi . Down- 
stream diffusion yields the dominant contribution 
to Tacc , with upstream particles requiring only 
small defiections (accomplished in short times: 
e.g., Illp from the shock normal in order to re- 
turn downstream. This automatically implies a 
hard lower bound to Tacc as Fi — > oo , since 
the downstream speed saturates at c/3. Effec- 
tively, particles can never be accelerated at rates 
much faster than their gyrofrequency. The limit 
translates to a comparable limit in the upstream 
fluid frame, which is often the observer's reference 
perspective, for example the interstellar medium 
surrounding a jet. This property follows from 
the connection between Lorentz transformations 
of times and energies, with the proper time of the 
particle being a Lorentz invariant. 

Hence, models of acceleration at relativistic 
shocks do not incur any increases to the energiza- 
tion rate other than the enhancement of the field 
by a single Lorentz boost. Maximum energies are 
then only explicitly weakly dependent on shock 
speeds. This implies that sites for cosmic ray ac- 
celeration generally must invoke higher environ- 
mental magnetic fields to effect higher maximum 
particle energies. 
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4. IMPACT ON UHECR PARADIGMS 

As an illustrative and topical case, the focus 
here is on the scenario that gamma-ray bursts 
are sites for the generation of ultra-high energy 
cosmic rays. The discussion will emphasize accel- 
eration issues, as opposed to source population 
considerations. 

4.1. Gamma- Ray Burst Applications 

The paradigm that gamma-ray bursts (GRBs) 
are responsible for UHECR production IfidlfiSlfiHj , 
has been quite topical over the last decade. 
Bursts are sufficiently energetic to amply satisfy 
cosmic ray energy budgets if their space density 
is not too sparse. Consequently the redshift dis- 
tribution of GRBs, not well-known at present but 
soon to be refined by the SWIFT mission, is an 
important constraint '61' on the ability of bursts 
to serve as sites for UHECR generation. Given 
the level of community interest, it is salient to 
assess the aforementioned acceleration results in 
the context of gamma-ray bursts. 

First, the discussions above indicate that the 
maximum energy of cosmic rays from bursts is not 
increased by subtle relativistic effects, and can 
be approximately estimated using standard non- 
relativistic shock theory with modification as per 
Eq. Q. It can be quickly deduced, by comparing 
the inverse gyrofrequency with subsecond burst 
variability timescales, that UHECRs can be gen- 
erated in GRBs if the magnetic fields inside the 
burst are of the order of i? ^ 10^ — 10''G inside the 
emission region. This is not dissimilar from field 
estimates (e.g., [SSEl) obtained by synchrotron 
radiation modeling of their prompt gamma-ray 
emission in the MeV band, so approximate con- 
sistency is achieved. It is not yet understood 
how such large fields arise in GRB shocks, though 
ideas of field amplification at shocks j51l48| have 
recently become prominent. 

The non-universality of the power-law index 
and its sensitivity to obliquity and the anisotropy 
of turbulent transport immediately indicate that 
GRB spectra should possess diverse indices. This 
should be manifested in the energy range above 
the 100 keV - 1 MeV peak of emission, and is in 
fact so in data taken from the EGRET experi- 



ment on the Compton Gamma-Ray Observatory 
(CGRO), where the half dozen or so bursts seen 
at high energies have a broad range of spectral in- 
dices 21 , namely a ~ 2 — 3.7 for dn/de^ cx £~" . 
This result suffers from limited statistics due to 
(i) the nature of bursts, and (ii) to EGRET's field 
of view being more limited than that for BATSE, 
the principal GRB experiment on CGRO. The 
upcoming GLAST mission will provide a more 
refined determination of the distribution of burst 
indices above 30 MeV after its launch in 2007. 

At the same time, in order to emit the radiation 
detected, GRBs must have underlying electron 
distributions that are relatively flat. By exten- 
sion, since the electron and cosmic ray ion dis- 
tributions trace each other in theories of shock 
acceleration where the diffusive mean free path 
is dependent only on rigidity, then the ion index 
must lie in the range cr ~ 3 — 6 . This is a spec- 
tral constraint commensurate with that imposed 
by the observed E'^ UHECR distribution [PISS] . 
As the UHECR spectrum results from a convolu- 
tion of a host of sources, modulo unknown prop- 
agation effects, one expects that the flatter dis- 
tributions will dominate, so there is a satisfying 
consistency between UHECR ion and GRB pho- 
ton spectra. Since GRB shocks are believed to 
be ultrarelativistic, the acceleration results dis- 
cussed above (as exemplified in Fig. [SJ indicate 
that strong cross field diffusion (i.e., k±_ ^ K|| ) 
will be necessary in the majority of bursts, if their 
shocks are oblique, as is highly likely. This prop- 
erty is required to provide a cosmic ray distribu- 
tion at least as flat as the observations, and in- 
deed efficiently generate cosmic rays in sufficient 
numbers. Note that the same would apply to jets 
in active galaxies such as blazars, as an alterna- 
tive source of UHECRs. 

An acceleration issue not addressed above con- 
cerns the shape of the particle distributions at 
thermal and slightly suprathermal energies. This 
is essentially an injection or dissipational heat- 
ing issue that is readily probed for electrons by 
the spectrum of prompt GRB emission. Bar- 
ing & Braby pursued a program of spectral 
fitting of GRB emission using a linear combina- 
tion of thermal and non-thermal electron popu- 
lations. These fits demanded that the prepon- 
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derance of electrons that are responsible for the 
prompt emission reside in an intrinsically non- 
thermal population, strongly contrasting particle 
distributions obtained from acceleration simula- 
tions. This result implies a conflict for acceler- 
ation scenarios where the non-thermal electrons 
are drawn directly from a thermal gas (the virtu- 
ally ubiquitous case), unless radiative efficiencies 
only become significant at highly superthermal 
energies. Another potential caveat is that strong 
radiative self- absorption could be acting, in which 
case the GRB spectral probe is not sampling the 
thermal electrons. Considerable work is needed 
to resolve this issue to ascertain whether an ac- 
celeration paradigm can be truly consistent with 
the GRB emission that is seen. 

5. OUTLOOK 

The UHECR field is clearly anticipating the 
next generation of observational data from the 
Auger experiment, just around the corner. In the 
meantime, theorists will continue to develop their 
models and hone their understanding. In terms 
of acceleration theory, at least two key develop- 
ments can be expected in the coming years. 

First, with the advent of rapidly increasing 
computational capability, reliance on analytic 
and Monte Carlo techniques is no longer essential, 
and these acceleration approaches will be sup- 
planted in part by modeling by fully 3D plasma 
codes, namely particle-in-cell (PIC) simulations. 
Such simulations have generated interesting re- 
sults in the last decade, but have been hampered 
|4()l42j by restricted dimensionality imposed by 
CPU memory and speed limitations. In the last 
2-3 years, results from 3D codes have achieved in- 
creasing visibility, particularly for modeling rela- 
tivistic shocks |62I5()| . Their key impediment is 
having to treat widely disparate inertial scales as- 
sociated with the proton to electron mass ratio. 
Their virtue is the accurate modeling of parti- 
cle heating, coherent acceleration effects and field 
amplification in the shock layer. One can expect 
many interesting developments and results from 
PIC codes in the coming decade. 

Monte Carlo techniques will still continue to be 
powerful tools, due to their capability of handling 



large dynamic ranges in spatial and momentum 
scales. Thus they are ideally suited to nonlin- 
ear acceleration scenarios, which can impact the 
interpretation of UHECR ion composition stud- 
ies, since nonlinear non-relativistic systems are 
well-known to preferentially accelerate ions with 
higher mass to charge ratios The applica- 

tion of the Monte Carlo approach to nonlinear, 
relativistic shocks has so far been very limited 
|28| . and this is territory ripe for investigation. It 
is anticipated that nonlinear effects will be more 
subtle for the ui ~ c domain, since particle dis- 
tributions are sensitive to the shock speed, obliq- 
uity and the type of scattering that operates. In 
addition, the shock hydrodynamics are dependent 
on the plasma anisotropy, and this is inextricably 
connected to ui , Oeni and 6'scatt also. Hence, 
research in the near future should focus on elu- 
cidating the interplay between these ingredients. 
Such simulations also need to explore ways to ad- 
dress injection and heating of thermal species in 
a more consistent manner, including the influence 
of electric potentials in the shock layer. 
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